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We calculate two-pion correlation functions at finite rapidities based on a hydrodynamical model which does 
not assume explicit boost invariance along the collision axis. Extracting the HBT radii through y} fits in both 
Cartesian and Yano-Koonin-Podgoretskii parametrizations, we compare them with the experimental results ob- 
tained by the PHOBOS. Based on the results, we discuss longitudinal expansion dynamics. 
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I. INTRODUCTION 

"Perfect fluidity" of the created matter at the Relativistic 
Heavy Ion Collider (RHIC) in BNL is one of the most excit- 
ing news in the field of high energy nuclear physics yj]. Ex- 
perimental results and their comparison with theoretical cal- 
culation reveal that the matter created at Au+Au collisions 
should be something like a liquid of quarks and gluons, un- 
like a gas of almost free partons as naively expected |2|]. One 
of the strong evidences of this finding is an observation of 
large elliptic flow (V2) and its agreement with a perfect fluid- 
dynamical calculation @]. In order to reproduce the experi- 
mental result with such models, an equation of state assuming 
partonic state at high temperature and a phase transition and 
rapid thermalization time (To < lfm/c) are required yfl. Now 
the hydrodynamic model based on numerical solutions of the 
relativistic hydrodynamic equation for perfect fluid becomes 
an indispensable tool for theoretical analyses of the relativistic 
heavy ion collisions. Furthermore, the model itself has been 
becoming more sophisticated to reproduce new experimental 
data with high statistics. Current most sophisticated ones are 
full three-dimensional (solving hydrodynamic equation with- 
out any symmetry) hydrodynamic expansion followed by a 
hadronic cascade Jj,|5|]. These models can reproduce most of 
soft hadronic observables. Especially, simultaneous descrip- 
tion of particle ratios, transverse momentum spectra and ellip- 
tic flow is possible with such hybrid models. 

However, there are still some insufficient ingredients in the 
hydrodynamic analyses. First, we don't have a reasonable ini- 
tial condition derived from the first principle. Recently, the 
Color Glass Condensate (CGC) has been proposed as a suit- 
able initial condition for relativistic heavy ion collisions JfJ. 
This picture has been examined as an initial condition for a 
hydrodynamic model in Ref. [7] and found to give a good 
description of some observables in the case of fully hydro- 
dynamic description of the collision process. However, this 
initial condition fails if one takes hadronic dissipation into ac- 
count [4]. This fact suggests there is still an open space for 
dissipative partonic phase, or improvement of the initial con- 
dition. 

Second, the equation of state (EoS) of the QCD matter has 
not been fully understood yet. Since one of the most impor- 
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tant merit of using hydrodynamic models is that it can be di- 
rectly related to the EoS, detailed information on the EoS for 
all region of temperature and baryonic chemical potential is 
indispensable. As for the RHIC energy, net baryon number 
observed at midrapidity is small enough to neglect it [8]. Nev- 
ertheless, EoS at finite baryonic chemical potential may play 
important role at forward rapidity region and heavy ion col- 
lisions at lower energies. Because of a well-known difficulty 
of lattice QCD at finite baryonic chemical potential [9], lattice 
QCD calculation has not provided the complete solution yet. 
For vanishing baryonic chemical potential, the lattice equa- 
tion of state clearly shows a different behavior from the free 
parton gas [ 10] and a lattice-inspired EoS is implemented in a 
hydrodynamic calculation ll ill . 

At last, in spite of the success in most of soft observables, 
results of the two-pion momentum intensity correlation from 
such hydrodynamical models do not agree with experimen- 
tal data yet. According to the symmetry of the wave func- 
tion of two identical bosons, the two-particle correlation func- 
tion can be related to sizes of the source from which parti- 
cles are emitted. This fact is known as Hanbury Brown-Twiss 
(HBT) effect. Because it concerns with source sizes, which 
depend on momentum of particle pairs due to collective flow, 
the pion correlation function should be a diagnostic tool for 
the space-time evolution of the matter. Since the disagree- 
ment was firstly found with a (2+l)-dimensional model with 
boost invariance along the collision axis lfl2tl . man y ex ten- 
sions such as an explicit longitudinal expansion II 1 3L 1 1411. in - 
corporating chemical freeze-out 03 > chiral model EoS lfl5ll . 
opaque source 111611. fluctuating initial conditions and contin- 
uous freeze-out 111711 have been examined and the discrepancy 
has been getting improved, but the situation is still unsatisfac- 
tory. There are various possibilities for further improvements. 

So far discussion on the HBT radii at the RHIC has been 
limited to midrapidity because of acceptances of two experi- 
ment, STAR and PHENIX. PHOBOS also has measured the 
two-pion correlation function. By virtue of the wider accep- 
tance of the detector, measurements at non-zero rapidity win- 
dows can be done, and the data are now available in Ref. OH- 
For analyses of such data in terms of the Cartesian parame- 
terization O9II20I1 . it should be noted that there exists an ad- 
ditional HBT radius called "out-long cross term" [21] which 
vanishes at midrapidity due to the symmetry. This radius con- 
tains information on the correlation between freeze-out points 
on the transverse plane and those on the longitudinal direction. 
Hence, it is expected that this quantity is sensitive to longitu- 
dinal expansion dynamics beyond boost-invariant approxima- 
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lion. Similar consideration also holds for the Yano-Koonin- 
Podgoretskii parametrization which has three radius parame- 
ters and one velocity parameter called YK velocity fl22l l23ll . 
The PHOBOS data also provide rapidity dependence of the 
YKP radii and YK velocity 111811 . which may impose a restric- 
tion on the longitudinal expansion dynamics. Indeed, the ini- 
tial matter distribution as an input for hydrodynamic calcula- 
tions has not been fixed yet. This is indicated by Hirano in 
Ref. f24ll . in which two different initial energy density distri- 
butions can provide reasonable agreement with experimental 
data of pseudorapidity distribution of charged hadrons mea- 
sured in BOA GeV Au+Au collisions at RHIC. 

In this work, we employ two different initial energy density 
distributions for the hydrodynamic equations, as in Ref. l24tl . 
We focus our discussion on central collisions. Both of them 
are so tuned that they reproduce the pseudorapidity distribu- 
tion of charged hadrons measured in the most central events 
at 200A GeV Au+Au collisions. Then we compare the space- 
time evolution and shape of the freeze-out hypersurface of the 
fluids and see how the difference in the initial condition is re- 
flected onto them. We calculate the two-pion correlation func- 
tion as the most promising experimental observable to see the 
difference. Extracting the HBT radii through Gaussian fits, we 
compare them with the experimental results and discuss the 
transverse momentum and rapidity dependence of the HBT 
radii. In the next section, we briefly review the hydrodynam- 
ical model used in this work. Initial conditions are given in 
Sec|III] In Sec. lIVI we show numerical solutions of hydrody- 
namical equations for the initial conditions given in Sec. [Ill] 
Results for the HBT radii as compared with the experimental 
data are given in Sec.[V] Section|VT]is devoted to a summary. 



II. HYDRODYNAMICAL MODEL 

The basic equation of hydrodynamical models is the 
energy-momentum conservation law 



d M T^ = 0, 



(1) 



where is the energy-momentum tensor. For a perfect fluid, 

T" v = (e + P)u fJ u v -Pg fjV , (2) 

where g*™ = diag(+, — , — , — ) and e, P and u p are energy den- 
sity, pressure and the four velocities of the fluid, respectively. 
If one takes a conserved charge i such as baryon number and 
strangeness into account, the conservation law 







(3) 



is added to be solved. Providing an EoS P = P(£,n,), one can 
solve these coupled equations numerically. 

In this work, we consider the baryon number charge as a 
conserved charge and adopt an equation of state which ex- 
hibits a first order phase transition on the phase boundary in 
T — /jb plane from the free massless partonic gas with three 
flavors to the free resonance gas which consists of hadrons 
except for hyperons up to 2 GeV/c 2 of mass with excluded 



volume correction|49j]. See Ref. H25H for the detail. The criti- 
cal temperature T c at vanishing baryonic chemical potential is 
set to 160 MeV. This model of current use is basically same 
with the one used in Refs. lfl3l [l6ll . 

Putting z-axis as the collision axis, we use a cylindrical co- 
ordinate system (t,T| s , r, ())) as follows; 



t = Tcoshr| s , 
z = Tsinhr| s , 
r x = rcoscj), 
r y = rsin(|). 



(4) 
(5) 
(6) 
(7) 



Here, X = \Jt 2 —z 2 is the proper time and r| s = l/21n[(f + 
Z) I (t — z)] is the space-time rapidity. Since we focus on central 
collisions, we assume azimuthally symmetric system. Then, 
by virtue of u^if = 1 , the four velocities are given in terms of 
a longitudinal flow rapidity Fl and a transverse flow rapidity 
F T as 



u z = cos!i(Fl — T| s ) coshFj, 
ifl s = sinh(F L - r| s ) coshF T , 
u r = sinhFL. 



(8) 
(9) 
(10) 



To solve the equations numerically, we employed a method 
based on the Lagrangian hydrodynamics which traces flux of 
the current. The numerical procedure is described in Ref. ll26Tl . 
For a treatment of the first order phase transition, we intro- 
duce a fraction of the volume of the QGP phase to express the 
energy density and net baryon number density at the mixed 
phase [25]. In this algorithm, we explicitly solve the entropy 
and baryon number conservation law. We checked that these 
quantities are conserved throughout the numerical calculation 
within 1% of accuracy for a time step 8x = 0.01 fm/c, by 
choosing proper mesh sizes of r| s and r directions. 



III. INITIAL CONDITIONS 

Firstly, we choose an initial proper time as Xo = 1 fm/c. 
Initial values for other variables are given on this hyper- 
bola. Longitudinal flow rapidity is set to the Bjorken's scaling 
ansatz Fl = r| s i27il . Transverse flow is simply neglected at 
the initial proper time 15011 . For the matter distributions, we 
assume that the energy and baryon number density are pro- 
portional to the number of binary collisions. Hence, for the 
Woods-Saxon profile of nucleon density in nuclei 



Pw0,z) 



Po 



,(V^+7-*)/? + l' 



(11) 



where R = 1.12A 1 / 3 - 0.86A -1 / 3 fm is the radius of the nu- 
clear with mass number A, Z, = 0.54 fm is the surface diffuse- 
ness and po is the normal nuclear matter density, the density 
of binary collisions at vanishing impact parameter is given by 
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(12) 



Kenji Morita 



3 



with Go being the total inelastic nucleon-nucleon cross section 
which is absorbed into the proportionality constant between 
«bc and matter distributions. 

Then, the energy density distribution is parameterized with 
a "flat+Gaussian" form, 



TABLE I: Parameters in initial matter distributions 



e(xo,ii s ,r) =£oexp 



(|r| s |-Tl s o) ; 



8(hls|-THo) 



«BCM- 



(13) 

Here, «bc('") is the normalized density of binary collisions 
(fT2] >. £o the maximum energy density, and r| s o and are 
parameters which determine the length of the flat region and 
width of the Gaussian part, respectively. Similarly, the net 
baryon number density distribution is parameterized as 



n B (lQ^s,r) =«Bo|exp 
(rUo-rUo) 2 



(hs|-TlsD) 



21 



2<& 



exp 



6(r| s o - h s 



OChsl-Tlso) 

Mc(r), (14) 



where «bo is the maximum net baryon number density and 
T| s d and o S D are the shape parameters as in Eq. d!3t . 

To calculate final particle distribution, we use the Cooper- 
Frye prescription 12811 . The pseudorapidity distribution for a 
particle species i is given by 



dNi di 
~dr\ ~ (2k) 2 Jo 



' kM f 

dkt J k ■ dof(k ■ u,T,/lib), 



(15) 



where tf is the momentum of thermally produced particles i 
with di being the degree of freedom, pseudorapidity r| defined 
byri = l/21n[(|k|+fc z )/(|k|-fc z )] wdf(k-u,T,ftB) the equi- 
librium distribution functions. We take into account not only 
directly produced particles but also resonance decay contribu- 
tions. The freeze-out hypersurface E is chosen to be a con- 
stant temperature one, T = Tf = 140 MeV. Here, we assume 
thermal and chemical freeze-out occur simultaneously. Since 
experimental data of particle yields can be well described by 
the statistical model with high chemical freeze-out tempera- 
ture close to T c Il29ll . we cannot reproduce the correct particle 
yields with this lower freeze-out temperature. However, in hy- 
drodynamic analyses, thermal freeze-out temperature is sensi- 
tive to k t spectra which mainly reflects transverse expansion. 
In this calculation, we set the freeze-out temperature so that 
pion k t spectrum is roughly reproduced in IC B and set the 
same freeze-out temperature for IC A and IC B. Note that the 
freeze-out temperature depends on the choice of the transverse 
profile of the initial matter distribution because a steeper pres- 
sure gradient yields larger transverse flow. For example, even 
Tf ~ 150— 160 MeV is possible with an initialization based 
on pQCD+saturation model 13011 . Our value is only slightly 
different from Ref. ll24ll where the initial profile is very sim- 
ilar. In order to reproduce both of the particle yields and the 
k t spectra in dynamical regimes, one should introduce sepa- 
rate freeze-out temperatures ffl l3lll or go to hybrid approach 
iB HI HH [IH] . In this work, however, our main argument will 
not be so affected by the description of the freeze-out because 
we focus on longitudinal expansion. 
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In Table |U two sets of the initial parameters are listed. The 
corresponding initial energy density distributions, resultant 
pseudorapidity distributions and transverse momentum distri- 
butions are illustrated in Fig. Q] [2] and [3] respectively. We have 
chosen two initial conditions both of which reproduce the ex- 
perimental data of PHOBOS on pseudo-rapidity distribution 
of charged hadron ll34ll . PHENIX on transverse momentum 
distributions of n~, K and p lj35ll . and BRAHMS on rapidity 
distribution of net protons [36]. These two initial conditions 
are characterized by two parameters, r| s o and o s o- One has 
small T| s o and large a s o, which we denote initial condition A 
(IC.A). The other which we represent IC.B, has an opposite 
feature; r| s o is large and o s o is small. The initial energy densi- 
ties are both much larger than experimental estimations (~ 5 
GeV/fm 3 ) based on Bjorken's formula ]2[L but note that £0 in 
Table|T]is not an average energy density but maximum energy 
density, which strongly depends on the profile of initial matter 
distributions ll37ll . We calculated the pseudorapidity distribu- 
tions for not only these two initial conditions but also inter- 
mediate ones by varying r| s o from 1.0 to 3.0 and found that 
they can also reproduce the experimental data by adjusting 
other parameters appropriately. Perhaps the best fit will exist 
in the middle of this parameter range HU. Here, we choose 
the extreme cases in order to see differences in the space-time 
evolution of the fluids originating from the difference in the 
initial conditions. 



IV. SPACE-TIME EVOLUTION OF THE FLUIDS 

Figures |4] and [5] show the space-time evolution of the tem- 
perature distributions and deviation from the scaling solution 
Yl = f\s> as a function of r| s at r = for various T, respec- 
tively. From these figures, we find that the space-time evo- 
lution at forward rapidity is quite different between IC.A and 
IC.B in spite of the fact that both solutions give similar pseu- 
dorapidity distributions of hadrons. In IC.B, sharp decrease 
of temperature which is identical to steep pressure gradient at 
the forward rapidity causes rapid acceleration of the longitu- 
dinal flow at the edge of the fluid. On the other hand, in IC.A, 
pressure gradient is rather gradual. Hence, resultant deviation 
from the scaling solution is smaller. Because pressure gradi- 
ent exist at smaller T| s in IC.A, however, such deviations take 
place at r| s ~ 1 while the flow keeps the scaling solutions up 
to r| s ~ 2 in IC.B. This fact explains slightly larger 80 in IC.A 
since faster longitudinal expansion than the scaling expansion 
pushes entropy per unit rapidity to forward rapidity lfl3l \3% . 

Although Figs. |4]and|5]show that there exist differences be- 
tween IC.A and IC.B in the space-time evolution, it is not triv- 
ial that such differences can survive at the freeze-out hypersur- 
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FIG. 1 : Initial energy density distributions 



FIG. 4: Space-time evolution of temperature distributions at r — 0. 



800 



z 

T3 




-2 2 

Pseudorapidity r\ 

FIG. 2: Pseudorapidity distribution. The solid line and dashed line 
stand for the initial condition A and B, respectively. Experimental 
data measured by PHOBOS are taken from Ref.| 34] 
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FIG. 3: Transverse momentum distribution of identified negatively 
charged hadrons. Spectra for kaons and anti-protons are scaled by 
factor 0.1 and 0.01 for a clear comparison of the slopes, respectively. 
Experimental data measured by PHENIX are taken from Ref. 13511 . 
Identification of symbols is same as Fig. [2] 



faces. Since hadrons strongly interact and have information 
only at the thermal freeze-out, differences on the freeze-out 
hypersurfaces are necessary to find the signature in hadronic 
experimental observables. 

We show freeze-out proper time Tf of the all fluid elements 
in Fig. [6] This characterizes the shape of the freeze-out hyper- 
surface, which is expected to affect the HBT radii. In Fig. [6] 
we can see that the system expands in the transverse direction 
in both of the fluids. Due to the same transverse profile, there 
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FIG. 5: Deviation from the scaling solution at r = 0. 
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FIG. 6: Freeze-out hypersurface of the fluids. 
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FIG. 7: Deviation from the scaling solution on the freeze-out hyper- 
surfaces. 



is no apparent difference in the transverse direction. On the 
other hand, the shape of the hypersurface in the r| s direction 
shows some variations. In IC.B, expansion appears and the 
freeze-out proper time is mostly constant in the broad range 
of T| s while it moderately decreases with r| s in IC.A. This is a 
consequence of the different longitudinal flow profile (Fig. [5]). 
We also plot the deviation from the scaling solution at the 
freeze-out in Fig. [7] The large deviation seen at forward ra- 
pidity in IC.B (Pig. O survives at the freeze-out. We will see 
how these differences affect the HBT radii in the next section. 



V. HBT RADII 

A. Two-pion correlation function 

Assuming that the source is completely chaotic, we can cal- 
culate the two-particle correlation momentum intensity corre- 
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lation function through this formula |4C 

\Kq,K)\ 2 



C 2 (q,K) = 1 



7(0,^)7(0^2)' 



(16) 



where q = k\ — ki is the four-relative momentum and K — 
\/2(ki + kj) is the four-average momentum, with k[ being 
on-shell momentum of emitted pions. The interference term 
I(q,K) can be chosen as 

I(q,K) = j K-dae^fiu-KJ), (17) 



so that 7(0, A;,) reduces to the Cooper-Frye formula 14111 . 

Experimentally, the two-pion correlation function is de- 
fined as 



C 2 (q) - 



B(qY 



(18) 



where A(q) is the measured two-pion pair distribution with 
momentum difference q and B(q) is the background pair dis- 
tribution generated from mixed events. Momentum accep- 
tances are imposed separately in the numerator and the de- 
nominator. Accounting for the large acceptance in the PHO- 
BOS experiment, 0.4 < Y %% < 1 .3 for three Kj bins and 0. 1 < 
Kj < 1 .4 GeV/c for three rapidity bins, we integrate the cor- 
relation function as follows: 



ftf dY™\I(q,K)\ 2 
/ 1 i^7(0,fei)7(0,fe 2 )' 
J^dK T K T \I(q,K)\ 2 
Jo 1 ! 4 dK T #r/(0,*i)/(<U 2 )' 



C(q;K T ) = 1 



C(q-Y KK ) = l 



(19) 
(20) 



For simplicity, we consider only directly emitted pions and 
neglect resonance decay contributions. 



B. Kj dependence of the HBT radii in the Cartesian 
parametrization 

Physical meaning of the HBT radii depends on the choice 
of three independent components of the relative momentum 
q. The most standard choice is the so-called Cartesian Bertch- 
Pratt parametrization |[l9|,[2(J q = (<?out,<7side,<7iong) in which 
"long" means parallel to the collision axis, "side" perpendicu- 
lar to the transverse component of the average momentum Kx 
and "out" parallel to Kx- In the case of azimuthally symmet- 
ric system as considered here, one can put Kx = (7Tx,0) so 
that <7 ou t = q x and g s id e = q y - Note that <7i ong = q z . Then, the 
Gaussian form of the two-pion correlation function is given as 



C 2 fit(q) = l+Xexp(-q 2 om R 2 



n 1 d2 2 d2 

'/side /l side '/lonsrMona 



-2q out q\ on gR \) 



(21) 



The HBT radii 7?, can be extracted by a % 2 -fit to the above 
fitting function. For a chaotic source, the chaoticity param- 
eter X should become unity. However, experimentally ob- 
served chaoticity is smaller than 1 because of such contribu- 
tions as long-lived resonance decay II42I1 . Here we fix X = 1 



in the Gaussian fit to the calculated correlation functions with 
Eqs. <|T9j and ©. 

By expanding the correlation function (fTol l for q ■ x <C 1 , 
the size parameters 7?, can be related to second order moments 
of the source function 02111 . In the Cartesian parametrization, 
taking the longitudinal co-moving system (LCMS) makes the 
expression simple; 



K ut = ((f x -Vi_f) 



r\2\ 



= (r\ 2 ) 

^side = (f> 2 ) 



■2p ± (fJ) + Pi(f 2 ). 



(22) 
(23) 
(24) 
(25) 



where 



(AW) 



J L k-daf(u' 



kJ)A{x) 



J L k-daf(u- 



kj) 



(26) 



x = x-{x), and P ± = £ T /7i k . Hence, 7? ou t. #side and 7?i ong 
can be interpreted as a mixture of thickness of the source 
and emission duration, transverse source size and longitudi- 
nal source, seen from the LCMS, respectively. Validity of 
these expressions for a hydrodynamical model is discussed 
in Ref . B43I1 . Although they have been shown to be good ap- 
proximations, it is also pointed out that there are still some 
discrepancies and one should use fitted HBT radii for com- 
parison with the experimental data which are obtained from 
the fit H. 

Figure [8] shows results for the four HBT radii compared 
with the experimental data measured by PHOBOS B18I1 . For 
comparison of the initial conditions, any qualitative and quan- 
titative difference cannot be seen in R out and R s id e , as expected 
from Fig. [6] 7?i ong of IC.A is about 1 fm smaller than that of 
IC.B. This can be considered as a consequence of the fact that 
the deviation from the scaling solution at small r| s is larger in 
IC.A, because faster flow causes more thermal suppression of 
the emission region fl3ll . For these three radii, our calcula- 
tion cannot reproduce the experimental results and show sim- 
ilar behavior with other perfect fluid dynamical calculations 
of Ref.lH [H [3 [H]. Especially R long shows the largest 
deviation from the experimental data, although calculation is 
improved by the longitudinal expansion without explicit boost 
invariance 11311 1411 . In the bottom of Fig. [8] result of the out- 
long cross term is presented. Reflecting the uniform shape of 
the freeze-out hypersurface in Fig. [6] the value of R^ of IC.B 
is smaller than that of IC.A. At the lowest Kj bin, the differ- 
ence is about 4 fm 2 . Unfortunately, experimental uncertainty 
is still too large to distinguish which initial condition is fa- 
vored. However, it should be noted that both of two results 
agree with the experimental data, in spite of the disagreement 
of other radii. 
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FIG. 8: Kj dependence of Cartesian HBT radii. Closed squares and 
open squares denote our results for IC.A and IC.B, respectively. Ex- 
perimental data are taken from Ref. fljjll . Error-bars for the experi- 
mental data are statistical only. 



C. Rapidity dependence of the HBT radii in the YKP 
parametrization 

In the YKP parametrization, three independent components 
of the relative momentum q are q± = Jq^ + q%, q\\ = q z = 
<7iong and q x = E\ — E^. Then, the Gaussian fitting correlation 



function is given as 

C 2Y Kp(q) = 1 + ^exp \-R\q\ - R\(q\ - q\) 
-(Rl+Rl)(q-Uf 



(27) 



where V = y( 1,0,0, vyk), Y= 1/y 1 ~ v yk and v 'yk is the 
fourth fitting parameter called YK velocity. The three HBT 
radii, R±, and R T are invariant under a longitudinal boost. 
Physical meaning of the parameters can be given in a similar 
manner [23] and becomes the simplest as follows, if one adopt 
the YK frame where vyk = 0, 



Rj 



{r~y 2 )=R 



(z 2 ) 



(?) 



long' 



(28) 
(29) 
(30) 



The main advantage of using YKP parametrization is that the 
three HBT radii directly give the transverse, longitudinal and 
temporal source size, that are seen from the YK frame. How- 
ever, one should note that the latter two, (l29l and (f30b , are 
approximate expressions which hold only if the source is not 
opaque 114311 . Hence, and R x cannot be always regarded 
as the source sizes in the presence of strong transverse flow 
which makes the source highly opaque lfl6ll . The general ex- 
pression of vyk is complicated one [23] but it can be regarded 
as a longitudinal flow velocity of the source measured in an 
observer's frame. 

We plot results of HBT radii for the YKP parametriza- 
tion in Fig. [9] Though PHOBOS measures only at small 
values of rapidity, we calculate the HBT radii for Y m = 
0.602,0.877, 1.122, 1.5, 2.0,2.5, 3.0,3.5 and 4.0 and show the 
results as a prediction. For comparison between IC.A and 
IC.B, R± seems to barely reflect the uniform structure along 
T| s direction in IC.B. While R^ shows the difference of order 1 
fm at small rapidity coming from the deviation of the scaling 
solution as well as in the third panel of Fig.[S] R T shows little 
difference but agrees with the experiment. Large experimental 
errors will be due to the known difficulty of the limited kine- 
matic region in the YKP parameterization ll45ll . Because of 
the large Kj window of the data, it is difficult to estimate geo- 
metrical opacity effect on R T . If we assume this effect is small, 
a possible origin of the deviation of our result from the data is 
larger emission duration. Some model calculations based on 
source parametrization [46] and parametric exact solution of 
hydrodynamics [47] show very small emission duration time, 
0-2 fm/c in agreement with data on the bottom panel of Fig. [9] 
We cannot see any significant differences in the HBT radii at 
forward rapidity expected from Figs. [6] and [7] which display 
the differences of the source shape and the longitudinal flow. 
This will come from the fact that the number of produced par- 
ticles is larger at late freeze-out proper time in the case of the 
current freeze-out condition j43ll . 

Finally the Yano-Koonin rapidity Fyk = l/21n[(l + 
vyk) /(l — vyk)] is shown as a function of Y nK in Fig. [TO] Both 
of results from IC.A and IC.B surprisingly agree with the ex- 
perimental data and show no difference between the two. At 
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FIG. 9: HBT radii for the YKP parametrization. The identification 
of the symbols is the same as in Fig. [8] 



forward rapidity region, our results show deviation from the 
infinite boost invariant case, which is indicated by the straight 
line. Although our solutions of longitudinal flow show devia- 
tion from the scaling solution (Figs.|5]and|7]i, the result would 
have to larger Fyk than a given Y nn if Fyk correctly represents 
the longitudinal source velocity. Hence, this deviation will be 
caused by the finite size effect [43] which becomes more sig- 
nificant at forward rapidity rather than the difference in the 
flow velocity. 



VI. SUMMARY 

In summary, we calculated the two-pion correlation func- 
tion for two sources which are given by a hydrodynamical 
model without explicitly boost invariance along the collision 
axis. The two initial conditions are so chosen that both of 



a, 

> 2 



IC.A 
IC.B 
PHOBOS n 
PHOBOS 7t + 



FIG. 10: The Yano-Koonin rapidity Fyk- The identification of the 
symbols is the same as in Figs.[8]and[9] The solid line indicates the 
case of the infinite boost-invariant source. 



them give consistent pseudorapidity distribution with the ex- 
perimental data and have different shape in the longitudinal di- 
rection. Other model ingredient, initial transverse profile from 
the binary collision model, scaling solution for initial longitu- 
dinal flow, vanishing initial transverse flow, EoS with first or- 
der phase transition and Cooper-Frye freeze-out prescription 
with 7f = 140 MeV are the same in the two solutions. We find 
that there exist some differences in the space-time evolution 
of the fluids in spite of the fact that both fluids give similar 
particle distribution. The HBT radii are extracted from the 
two-pion correlation functions and compared with the exper- 
iment. In the Cartesian parametrization, the out-long cross 
term which arises at nonzero rapidities shows a difference be- 
tween two initial conditions and the good agreements with the 
experimental data. The correlation function is also analyzed 
with the YKP parametrization. We find a small difference be- 
tween the two initial conditions in which reflects deviation 
from the scaling solution in the longitudinal expansion as well 
as /?iong in the Cartesian parametrization. Possible sources 
of this disagreement are followings: EoS of current use ex- 
hibits first order phase transtion which makes the lifetime of 
the fluid longer, and assumes hadronic states is in fully chem- 
ical equilibrium. It is known that both crossover EoS IU5I1 and 
incorporating chemical freeze-out 1 14] improve the lifetime 
then ^iong and R\\. We used the conventional Cooper-Frye 
prescription for the freeze-out. Improvement of the freeze- 
out prescription by continuous freeze-out ifTTll and hybrid ap- 
proach ]4j, |5[] can yield larger R^g but this may lead to larger 
R out because of extended emission duration. Nevertheless, as 
a transport calculation [48] shows, positive x — t correlation in 
the source function may resolve this problem. Finally, in spite 
of the disagreement of the HBT radii, the YK rapidity shows 
a good agreement with the experimental data. Our calculation 
predicts some deviations at larger rapidities from the infinite 
boost-invariant case. Hence, measurements at this region is 
needed for further understanding of the expansion dynamics. 
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